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Abstract 

Atomistic (atom-scale) and coarse-grained (meso-scale) simulations of 
structure and dynamics of poly-isoprene melts are compared. The lo- 
cal structure and chain packing is mainly determined by the atomistic 
details of the polymer architecture. The large-time dynamics encoun- 
tered in NMR experiments can be explained by meso-scale simulations 
including stiffness. The connecting link between the two scales is the 
stiffness which, although being a local property, influences strongly 
even the long-timescale dynamics. The standard reptation scenario 
fails to explain the observed dynamics. We propose strong reptation 
as a modified reptation scenario in which the local Rouse motion is 
absent. 
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1 Introduction 

The large abundance of polymers and their variety of applications 
make them an interesting target of study in theoretical material sci- 
ence. The understanding of the differences and similarities of the 
various materials is an important prerequisite for the goal of CAMD 
- computer aided materials design. The dream of chemical engineers 
would be to develop in the workstation the perfect material for a given 
purpose. 
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To come closer to this distant goal, much effort in various fields 
is necessary. As polymeric materials are characterized by the im- 
portance of various length scales, the understanding of the interplay 
between these scales is of utmost importance. Methods adapted to all 
relevant scales are needed from the experimental as well as from the 
theoretical or simulational viewpoint. In simulations, much work has 
been directed recently of the issue of coarse-graining, the mapping of 
simulations on different scales, in order to get a unified view of the 
arising scales [Q-[|]. In the present contribution we show that atom- 
istic simulations of irans-polyisoprene (PI) || can be mapped onto a 
simple bead-spring model incorporating exclude volume, connectivity, 
and stiffness 0, ^] . With this simple model we can investigate the 
long-time dynamics of chains in the melt. There is evidence for rep- 
tation which qualitatively changes with stiffness |^,|l(|. This connects 
to results of modern NMR experiments [jnj] which we can reproduce 
to a satisfactory extent (§|,[l^]. Thus, simulations on both scales and 
especially their connection reveal different important aspects of the 
system under study. 

The remainder of this contribution is organized as follows. In sec- 
tion |2] a short review about our recent results of atomistic simulations 
of Trans-PI is given. In section || results of simulations on the meso- 
scale level, where the polymer identity is put into a simple stiffness 
parameter, are presented. In the concluding section the important 
concept of stiffness is discussed and it is shown that stiffness is suffi- 
cient to allow a mapping for the polymer presented here. 

2 Local Structure and Dynamics - The 
atomistic scale 

The atomistic structure of oligomers of trans-1, 4-polyisoprene (cf. Fig- 
ure |l|) was investigated. For details of the simulations and the interac- 
tion potentials see refs. ||,||. Here we only note that the simulation box 
contains 100 oligomers of average length 10 monomers of pure trans- 
poly-isoprene which were pre-equilibrated using end-bridging Monte 



Carlo [13,14] at the ambient condition of 300 K or the elevated tem- 



perature of 413 K. All atomistic simulations are run at 101.3 kPa. 

The main result is, that the local structure is very important for 
properties on the sub-monomer length scale. Figure shows the local 
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mutual packing of chains by means of radial distribution functions of 
the different atoms in the melt. The carbons show a distinct peak at 
next neighbor contact. The hydrogens, in contrast, show very little 
structure. Thus, we encounter a carbon backbone with a surrounding 
"hydrogen cloud" . This can be taken as a first hint that not all details 
are necessary for every simulation. Still, the hydrogens had to be taken 
into account for the C— H vector reorientation. This is important for 
direct comparison to NMR experiments (below). The comparison to 
experimental raw data is always a good and necessary validation of 
simulation models [12|. 



The local packing is additionally reflected by the directionality of 
chains at contact. This is shown in Figure |3| by means of the spatial 
orientation correlation function of the double bonds 

OCF{r) := P 2 (r) = ^- 3cos 2 a(r) - 1 ^, (1) 

where a(r) is the angle between tangent unit vectors on two different 
chains. The distance r is measured between their centers of mass. 
The unit vectors may be defined in different manners, e.g. the double 
bonds in Figure [|| are one possibility to denote the direction of a 
monomer. 

Direct comparisons to experiments |l5|-[T7[| and simulations [18,19| 
on cis-polyisoprene proof that our model is realistic 



The reorientation of the hydrogens is subject to a two-stage pro- 
cess (cf. Fig. H|). The first stage is a fast initial decay on the time- 
scale of the segmental motion, i.e. the motion on the monomer or 
sub-monomer scale. On the second time-scale the reorientation of the 
whole oligomer becomes important. Overall, reorientation is moni- 
tored by the following correlation function 



3cos 2 $(t) - 1 ). (2) 



c 

^reor — \ ^ 

Here, is the angle by which a given bond vector reoriented in 
time t. The second Legendre polynomial is chosen as it is the rele- 
vant quantity in NMR measurements. The effectivity of the two stages 
may be measured by fitting double exponential decays to the obtained 
correlation functions. In doing so, satisfactory ageement with exper- 
iments could be achieved (cf. Table |T|) if we keep in mind that the 
investigated systems are not completely identical. The experiments 
were performed on mixtures with a high cis-PI content and longer 
chains. 
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Figure 2: Radial distribution functions of atomistic data of trans- 
polysoprene. a) Carbon— Carbon, b) Carbon— Hydrogen and c) 
Hydrogen— Hydrogen RDFs. The carbons show a stronger pronounced struc- 
ture. A correlation hole can be seen in subfigure c). 
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Figure 3: Atomistic mutual orientation of double-bonds (blue) in trans- 
polyisoprene in comparison to the fully flexible simple model of section ^ 
(red). The black line is for the atomistic vectors connecting C5 with Ci of 
the next monomer. 
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Figure 4: Reorientation correlation functions of the methylene C— H vectors 
at 413 K. Blue: C 5 -H, Black: Ci-H. 
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Table 1: Comparison of the experimental (cis-PI) and simulation (cis and 
trans-Pi) data for the efficiency of the initial stage of the reorientation pro- 
cess. The parameter a denotes to which value the reorientation correlation 
function (Eq. |^) decays in the short-time process, i.e. before the long-term 
exponential relaxation sets in. a lps is the value of 1 — C reor at lps. In the 
analysis of the simulations for cis-polyisoprene a stretched exponential sec- 
ond process was assumed. The experiments used a range of temperature 
between 283 K and 363 K [[DJ. The trans simulations were at 300 K || and 



the cis simulations at 363 K |19] 



In section 3.2 comparisons of the long-time dynamics of the corre- 
sponding simple model with experiments will be presented. 



3 Global Structure and Dynamics - Meso- 
scale simulations 

Simple polymer models allow one to investigate large-scale phenomena 
both in time and space which are not accessible by atomistic simula- 
tions. Therefore many researchers employed such models to look for 
rather generic polymer properties or dynamical concepts [pC|-p5| . One 
of these important concepts is reptation |2q, 27|. The reptation con- 
cept explains much of the molecular weight dependence of viscosity 



and the elastic and loss moduli [27]. However, this model originates 
in the simple Rouse model |2q] where no stiffness is included or it is 
subsumed into simple Kuhn blobs [29j. 

In order to look for the influence of stiffness, a three-body potential 
for stiffening the backbone has therefore to be introduced P, p|,|30|] 
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Figure 5: Static structure factor of model chains of length 50 monomers 
depending on stiffness. Stiffness is given in terms of the persistence length 
(see text). "Flex" means that no stiffening potential is introduced; excluded 
volume leads to a persistence length of l p — 1. 

The force constant l p in this choice of units has the same numeri- 
cal value as the resulting intrinsic persistence length (below), and is 
therefore denoted as l p . 

3.1 Statics 

In polymer melts excluded volume is commonly assumed to be screened 
out. Thus, we can expect that polymer chains behave as random 
walks ^7|. This is true at least on large length scales. All local inter- 
actions only result in local chain stretching. This is seen easily in the 
single chain static structure factors of model chains |7]] (cf. Figure [5]). 
The random walk appears in the fractal dimension of dt = 2 on the 
length scale bigger than the persistence length. The fractal dimen- 
sion expresses itself in the slope of the structure factor. The chains 
with stiffness bend over to a weaker slope indicating a smaller frac- 
tal dimension at higher \k\. The persistence length l p originates from 
the Kratky-Porod worm- like chain picture |3l]| ; l p measures the decay 
length of bond correlations along the chain backbone. Thus, static 
properties of different models and chains are easily mapped onto each 
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Figure 6: Static mutual chain packing of chains of different stiffness. The 
stiffer chains order more strongly parallel. Long range isotropy is preserved 
(OCF -> 0,r -> oo). 



other by the simple blob concept where all local interactions are put 
into one single length scale, which then allows a renormalization onto 
flexible chains with the blobs acting as coarse-grained monomers. The 
size of the Kuhn blob Ik and the persistence length l p differ only by 
the constant factor of 2 although they arise from different concepts. 

Still, even in chains of complete flexibility, i.e. no three-body po- 
tential, there is a small but visible alignment between neighbouring 
chains, similar to that shown in Figs. |3| and |6| on the length scale 
of about 3 chain diameters [24]. This is evidence that the very sim- 
ple ansatz of model chains only incorporating connectivity comes to 
its limits as soon as many-body interactions come into play. As we 
have seen in section ^ local interactions affect strongly the mutual 
packing of chains. If chain stiffness is increased, chain order becomes 
stronger Q (cf. Fig. ^) without leading to an overall nematic order. 
This ordering is a strictly local phenomenon which is proven by the 
fact that there is no chain length dependence whatsoever 
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Figure 7: Algebraic reorientation of segments of the chain backbone depend- 
ing on chain stiffness (Chain length: 200 monomers). 



3.2 Dynamics 

For simple flexible bead-spring chains the reptation concept was sue- 

.11,11] 



However, NMR 



cessfully validated by several simulations |22 
experiments of polymer reorientation show that this model cannot de- 
scribe all polymers satisfactorily. For poly-(dimethylsiloxane), PDMS, 
which is known to be very flexible the results of double quantum 
NMR can be explained by the simple reptation model |3!| whereas for 
polybutadiene this simple explanation shows strong deficiencies [11]. 
For this purpose the reorientation behaviour of entangled melts de- 
pending on stiffness was investigated J|] . With increasing stiffness the 
reorientation slows down tremendously going hand in hand with a de- 



creasing entanglement length and shrinking tube diameter [1C|. The 
reorientation of backbone segments is algebraic on short time scales 
(cf. figure 0). This algebraic dependence C reor oc i~ 7 is the same as 
found in NMR experiments. The exponents of 7 = — 7 and 7 = — ^ 
as seen in figure |8| are both found. However, there is a qualitative 
change in dynamics with chain length as the entanglement length is 
crossed. The exponent of 7 = — \ is not found in unentangled chains. 

In diffusion also the entanglement length scale is found. Chains 
longer than l e diffuse clearly slower than predicted by simple Rouse 
motion (Figure |9|) . According to the Rouse model the overall diffusion 
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Figure 8: Comparison of simulated (chains of length 200 and l p = 5; lines) 
and experimental (symbols) [jO]] reorientation correlation functions. For the 
experiments time-temperature superposition is assumed. 
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Table 2: Entanglement monomer number N e depending on persistence length 



l p , determined by the change of diffusion behavior with chain length |10]. 
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Figure 9: Diffusion coefficient of chains depending on chain length and stiff- 
ness. In the given representation Rouse behavior would give a straight hori- 
zontal line. 



is expected to be D oc N^ 1 . This would correspond to a horizontal 
straight line in figure ^. Reptation leads to D oc N~ 2 which we find 
for longer chains. The crossover point can be taken as a definition for 
the entanglement length |jl0| (Table ||) . We observe that the entangle- 
ment length decreases with increasing persistence length. l e and l p are 
neither independent nor linearly dependent on each other. This makes 
it impossible to renormalize stiff chains onto the simple bead-spring 
model, as only one length scale could be scaled away. The influence of 
stiffness survives on very long length scales. The entanglement length 
l e describes the anisotropy of motion of a chain due to the temporary 
network of its uncrossable neighbors. A polymer chain has to move 
predominantly along its backbone as transversal motion is hindered 
by the neighbors leading to an effective tube (cf. Figure |10|). This 
is measured by the correlation of directions of the chain (static u(t)) 
with the motion of the monomer in time (dynamic v(t — to)). 



C 



(l[3u(r)v(t-t )]) (4) 



This function measures the correlation between the static direction 
of a chain segment at a given point in time with the direction of its 
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Figure 10: Correlations of monomer motion with the direction of the back- 
bone depending on stiffness. Solid lines: t = t , dashed lines t = t + r/2 



dynamic motion in the time ther after. Thus, it shows that chain seg- 
ments move in the beginning preferably along their contour. This is 
exactly what reptation is about. Increasing stiffness supports this ef- 
fect as the stiffness suppresses the transversal motion even further. 
On short time-scales where in the standard reptation picture isotropic 
motion is still possible the local stiffness disallows this motion. Thus, 
the chains have to reptate from the very beginning. This is illus- 
trated by the mean-squared displacements of inner monomers ((x 2 ) 
cf. Figure 11). According to the standard reptation picture one ex- 
pects the scenario we find here only for flexible chains. On short 
time-scales a Rouse motion [28 is found ((x 2 ) cx t ' 5 ), then the Rouse 
motion is constrained to the tube ((x 2 ) oc i 0,25 ). After the internal 
degrees of freedom are relaxed the chain as a whole moves in the tube 
((x 2 ) oc i ' 5 ) and finally the chains reach the free diffusion ((x 2 ) oc t). 
The other extreme we now see for chains with a persistence length 
of five monomers. The first two dynamical regimes are missing com- 
pletely, as the Rouse motion is no more the correct description of the 
polymer especially on the short scales (see Figure |l2| ). This regime 
we like to call strong reptation. Similar results have been obtained 
by Morse p3|-|35[ for chains of stronger stiffness, who introduced the 
terms loosely entangled and strongly entangled for the different sys- 
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Figure 11: Mean squared displacements of central monomers in a chain de- 
pending on chain stiffness. For clarity lines indicating the slopes correspond- 
ing to t ' 25 and t 0,5 are shown additionally. 

terns, respectively. 



Note the intersections of the different curves in Figure 11. This 
shows that there are regimes where stiffer chains diffuse even faster as 
the entanglement length and tube diameter come down tremendously. 
This can also be seen if we look at dynamical structure factors Jl( 
which can measure the tube diameter. 



If the Rouse model was applicable in all sub-figures of Figure |12 
all curves would fall on top of each other as the Rouse modes 

1 ^ firp(i + 1/2) \ - 
with R{ the position of bead i and N the monomer number would 



be the true eigenmodes of the system and the time in Figure |12| is 
rescaled accordingly. However, we see that with increasing stiffness 
and chain length the model shows deficiencies. Because of stiffness 
the high modes (local motion) start to be affected and for long chains 
the entanglements hinder the low modes (large scale) motion. As the 
entanglement length shrinks with stiffness the Rouse regime is eaten 
up between the entangled motion on the one side and the local stiffness 
on the other. 
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Figure 12: Normalized correlation functions of Rouse modes X p of stiff and 
flexible chains of various stiffnesses (a) 1-50; b) 3-50; c) 5-200; d) 1-350). The 
notation means l p — N with persistence length and chain monomer number. 
The Rouse model breaks down with increasing length and stiffness. 



15 



4 Conclusions: Stiffness - A decisive 
characteristics 

Two different polymer models were introduced. A bead-spring model 
with stiffness and a fully detailed atomistic model. Both are validated 
against different experiments. The connecting link between the mod- 
els is the backbone stiffness which survives from the very local scale 
to the global scale on which often only entanglements are expected to 



be important. In Figure 13 one sees the success of the mapping. The 
atomistic chains at 413 K and model chains (also of length 10) are 
compared (one monomer to one monomer mapping). The mapping is 
accomplished by rescaling (squared) lengths with the mean-squared 
end-to-end distance. Time-scales are fixed by the center-of-mass dif- 
fusion. Then the figure shows the comparative reorientation of lo- 
cal monomer-to-monomer vectors. Thus, this study opens one possi- 
ble route to polymer coarse graining: Simulate an atomistic melt (of 
short chains) in full detail, there you can investigate all the very local 
observables (sub-monomer to monomer-scale) like radial distribution 
functions, orientation correlations up to structure functions and re- 
orientation times. Also one has to determine the persistence length 
of the polymer on this scale. This is then taken as an input to the 
model on the next length scale, so that part of the polymer identity 
is preserved and a tremendous simulation speedup is possible at the 
same time. At this scale, long-time dynamic phenomena like reptation 
and large scale static structure like the overall Gaussian random walk 
distribution can be looked into. 

However, one has to be very careful with this mapping. There are 
polymers for which this route is too simple-minded, especially if the 
monomer is strongly anisotropic or has special or bulky side-groups. 
More elaborate methods have to be applied then |l],§. Our method 
can be applied to dense melts of simple hydrocarbon polymers and 
allows in this case a very strong speedup and a look on real large 
scale phenomena. With the other methods typically an intermediate 
scale between the two presented scales here is introduced as additional 
length scales can become important. Still, for dynamic issues it is not 
possible even on the largest scale to explain everything with simple 
chains only, at least stiffness has to be taken into account. 
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Figure 13: The reorientation of atomistic chains and simple chains (both of 
length 10 monomers) in comparison. The red (Ci) and the blue (C2) line are 
atomistic vectors connecting the indicated carbons on neighboring monomers. 
The green line is the reorientation in the simple model with no added stiffness. 
The black line corresponds to the simple model with persistence length l p = 
1.5 similar to the persistence length of the polyisoprene model. To suppress 
end-effects the terminal monomers were not taken into account. 
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